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I. INTRODUCTION 

cn' 

The quantization of the gravitational interaction is a major open challenge in theoretical physics. This review 
. presents the status of the spin foam approach to the problem. Spin foam models are definitions of the path inte- 
' gral formulation of quantum general relativity and are expected to be the covariant counterpart of the background 
independent canonical quantization of general relativity known as loop quantum gravity jlHsl] . 
, ^ ■ This article concentrates on the definition of the recently introduced Engle-Pereira-Rovelli-Livine (EPRL) modeljl, 
Q and the closely related Freidel-Krasnov (FK) model [6] . An important original feature of the present paper is the 
explicit derivation of both the Riemannian and the Lorentzian models in terms of a notation that exhibits the close 
relationship between the two at the algebraic level that might signal a possible deeper relationship at the level of 
transition amplitudes. 

O ' We will take Plebanski's perspective where general relativity is formulated as a constrained BF theory (for a review 
^ ' introducing the new models from a bottom- up perspective see [3]; for an extended version of the present review 
?H including a wide collection of related work see 8]). For that reason it will be convenient to start this review by 
introducing the exact spin foam quantization of BF in the following section. In Section IIIII we present the EPRL 
model in both its Riemannian and Lorentzian versions. A unified treatment of the representation theory of the relevant 
—\ ' gauge groups is presented in that section. In Section HVl we introduce the FK model and discuss its relationship with 
I> ■ the EPRL model. In Section |V] we describe the structure of the boundary states of these model and emphasize the 
^ ' relationship with the kinematical Hilbert space of loop quantum gravity. In Section IVII we give a compendium of 
important issues (and associated references) that have been left out but which are important for future developpement. 
Finally, in section IVIII we present the resent encouraging results on the nature of the semiclassical limit of the new 
models. 
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II. SPINFOAM QUANTIZATION OF BF THEORY 



. 5^ r We will start by briefly reviewing the spin foam quantization of BF theory. This section will be the basic building 
' block for the construction of the models of quantum gravity that this article is about. The key idea is that the 
\^ ,. quantum transition amplitudes (computes in the path integral representation) of gravity can be obtained by suitably 
restricting the histories that are summed over in the spin foam representation of exactly solvable BF theory. We 
describe the nature of these constraints at the end of this section. 

Here one follow the perspective of . Let G be a compact group whose Lie algebra q has an invariant inner product 
here denoted () , and M a d-dimensional manifold. Classical BF theory is defined by the action 



S[B,uj] = y"(B AF 
M 



M), (1) 



where B is a g valued (d — 2)-form, w is a connection on a G principal bundle over JV[. The theory has no local 
excitations: all solutions of the equations of motion are locally related by gauge transformations. More precisely, the 
gauge symmetries of the action are the local G gauge transformations 

(5B ^ [B, a] , 5uj^ A^a, (2) 
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where a is a g-valued 0-form, and the 'topological' gauge transformation 

SB = dujri, 6io = 0, (3) 

where d^^ denotes the covariant exterior derivative and 77 is a g-valued O-form. The first invariance is manifest from 
the form of the action, while the second is a consequence of the Bianchi identity, di^F{u;) — 0. The gauge symmetries 
are so vast that all the solutions to the equations of motion are locally pure gauge. The theory has only global or 
topological degrees of freedom. 

For the moment we assume ^Vf to be a compact and orientable. The partition function, Z, is formally given by 

Z= jV[B]V[uj] exp{i I (BAF(w)}). (4) 

J J M 

Formally integrating over the B field in (|4]) we obtain 

Z = / <5(FH). (5) 



The partition function Z corresponds to the 'volume' of the space of flat connections on A4. 

In order to give a meaning to the formal expressions above, we replace the d-dimensional manifold A4 with an 
arbitrary cellular decomposition A. We also need the notion of the associated dual 2-complex of A denoted by A*. 
The dual 2-complex A* is a combinatorial object defined by a set of vertices v 6 A* (dual to d-cells in A) edges 
e e A* (dual to (d— l)-cells in A) and faces / £ A* (dual to (d— 2)-cells in A). In the case where A is a simplicial 
decomposition of 7W the structure of both A and A* is illustrated in Figures [U [21 and [3] in two, three, and four 
dimensions respelctively. 



FIG. 1: On the left: a triangulation and its dual in two dimensions. On the right: the dual two complex; faces (shaded 
polygene) are dual to 0-simpleces in 2d. 

For simplicity we concentrate in the case when A is a triangulation. The field B is associated with Lie algebra 
elements Bf assigned to faces / G A*. We can think of it as the integral of the (d— 2)-form B on the (d— 2)-cell dual 
to the face / e A*, namely 



Bf= J B. (6) 

(d-2)-ccn 

In other words Bf can be interpreted as the 'smearing' of the continuous (d— 2)-form B on the (d— 2)-cells in A. We 
use the one-to-one correspondence between faces / G A* and (d— 2)-cells in A to label the discretization of the B field 
Bf. The connection w is discretized by the assignment of group elements 6 G to edges e G A*. One can think of 
the group elements ge as the holonomy of uj along e G A*, namely 



9e 



Pexp(- / Lu), (7) 

J e 



where the symbol "P exp " denotes the path-order-exponential that reminds us of the relationship of the holonomy 
with the connection along the path e G A*. 
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FIG. 2: On the left: a triangulation and its dual in three dimensions. On the right: the dual two complex; faces (shaded wedge) 
are dual to 1-simpleces in 3d. 





FIG. 3: On the left: a triangulation and its dual in four dimensions. On the right: the dual two complex; faces (shaded wedge) 
are dual to triangles in 4d. The shaded triangle dual to the shaded face is exhibited. 



With all this the discretized version of the path integral (|4]) is 



Z(A)= / n dg, W dBfC^^f^f = / n '^Se n ^i9e,---gej 

'' eeA* /eA* '' eeA* /eA* 



(8) 



where Uf = ■ • ■ ge„ denotes the holonomy around faces, and the second equation is the result of the B integration: 
it can be thus regarded as the analog of ([5]). The integration measure dBf is the standard Lebesgue measure while the 
integration in the group variables is done in terms of the invariant measure in G (which is the unique Haar measure 
when G is compact). For given h ^ G and test function F{g) the invariance property reads as follows 

j dgF{g) = j dgF{g-') - J dgF{gh) - J dgF{hg) (9) 

The Peter- Weyl's theorem provides a useful formula or the Dirac delta distribution appearing in namely 

<5(.g)=^dpTr[p(g)], (10) 



where p are irreducible unitary representations of G. From the previous expression one obtains 

2(A) - Yl f li^de l[ dp, Tr [pfigl . . . )] . 

r.I „\^s f\'' »cA* /GA* 



(11) 



C:{p}-+{/} eeA* 



Integration over the connection can be performed as follows. In a triangulation A, the edges e G A* bound precisely 
d different faces; therefore, the (/e's in ([Tl]) appear in d different traces. The relevant formula is 



PtnviPir ■ ■ . Pd) / dge Pl{ge)(E)p2ige)(S■■■(E)Pdige)■ 
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For compact G it is easy to prove using the invariance (and normalization) of the the integration measure ([9]) that 
^inv — i^inv)^ ^^'^ projcctor onto Inv[pi ® P2 ® ■ • • ® Pd]- In this way the spin foam amphtudes of 50(4) BF theory 
reduce to 



ZBFiA) -En n 

C/:{/}-»P/ /eA* eSA* 



(13) 



In other words, the BF amphtude associated to a two complex A* is simply given by the sum over all possible 
assignments of irreducible representations of G to faces of the number obtained by the natural contraction of the 
network of projectors P[nv according to the pattern provided defined by the two-complex A*. 

There is a nice graphical representation of the partition function of BF theory that will be very useful for some 
calculations. On the one hand, using this graphical notation one can easily prove the discretization independence of 
the BF amplitudes. On the other hand this graphical notation will simplify the presentation of the new spin foam 
models of quantum gravity that will be considered in the following sections. This useful notation was introduced by 
Oeckl [l3, [HI and used in [l^] to give a general prove of the discretization independence of the BF partition function 
and the Turaev-Viro invariants for their definition on general cellular decompositions. 

Let us try to present this notation in more detail: The idea is to represent each representation matrix appearing in 
(fTTj) by a line (called a wire) labeled by an irreducible representation, and integrations on the group by a box (called a 
cable). The traces in equation pT|) imply that there is a wire, labelled by the representation pf, winding around each 
face / e A*. In addition, there is a cable (integration on the group) associated with each edge e G A*. As in (|13l) . 
there is a projector P^^^ is the projector into Inv[pi (g) p2 <X) • • • ® Pd] associated to each edge. This will be represented 
by a cable with d wires as shown in (jl4p . Such graphical representation allows for a simple diagrammatic expression 
of the BF quantum amplitudes. 



Pi P2 P3 



Pd 



,Pd) 



(14) 



The case of physical interest is d = 4. In such case edges are shared by four faces; each cable has now four wires. 
The cable wire diagram giving the BF amplitude is dictated by the combinatorics of the dual two complex A*. From 
Figure |3] one gets 



ZBFiA) = E n dp 




(15) 



The 10 wires corresponding to the 10 faces / e A* sharing a vertex v e A"^ are connected to the neighbouring vertices 
through the 5 cables (representing the projectors in (|13p and Figure [HI) associated to the 5 edges e S A* sharing the 
vertex v G A*. 



1. SU{2) X SU{2) BF theory: a starting point for 4d Riemanman gravity. 

We now present the BF quantum amplitudes in the case G = SU{2) x SU{2). This special case is of fundamental 
importance in the construction of the gravity models presented in the following sections. The product form of the 
structure group implies the simple relationship Zbf{SU{2) x SU{2)) = Zbf{SU{2)Y . Nevertheless, it is important 
for us to present this example in explicit form as it will provide the graphical notation that is needed to introduce the 
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gravity models in a simple manner. The spin foam representation of the BF partition function follows from expressing 
the projectors in (|15p in the orthonormal basis of intertwiners, i.e. invariant vectors in Inv[pi (X) • • • ® Pi]. From the 
product form of the structure group one has 



y y Ji jo Jj U jtjtjj jt 

TTTT llff WT 



Ji J2 J3 J4 Ji J2 J3 H 

\U/ NLy 



(16) 



where pf — jf ® J/^; ^-nd and are half integers labelling left and right representations of SU{2) that defined 
the irreducible unitary representations of G = SU{2) x SU{2), and we have used the expression of the right and left 
SU{2) projectors in a basis of intertwiners, namely 



il ^2 J3 ^4 Jl J2 J 3 J 4 

TTTT ' /TN 



(17) 



where the four-leg objects on the right hand side denote the invariant vectors spanning a basis of Inv[ji (8) ■ • • (8> J4] , 
and i is a half integer labelling those elements. Accordingly, when replacing the previous expression in (jlSp one gets 




(18) 



and equivalently 




(19) 



from which we finally obtain the spin foam representation of the SU (2) x SU (2) partition function as a product of 
two SU{2) amplitudes, namely 



6 




J4 'I 




(20) 



Extra remarks on four dimensional BF theory 



The state sum is generically divergent (due to the gauge freedom analo gou s to ([3])). A regularized version 
defined in terms of SUq{2) x SUq{2) was introduced by Crane and Yetter [Tsl. Il4l|. As in three dimensions, if an 
appropriate regularization of bubble divergences is provided, (1111) is topologically invariant and the spin foam path 
integral is discretization independent. 

As in the three dimensional case BF theory in any dimensions can be coupled to topological defects [l^. In the 
four dimensional case defects are string-like [IQl and can carry extra degrees of freedom such as topological Yang- 
Mills fields ■ The possibility that quantum gravity could be defined directly form these simple kind of topological 
theories has also been considered outside spin foams [ll] (for which the UV problem described in the introduction is 
absent) is attractive and should, in my view, be considered further. 

It is also possible to introduce one dimensional particles in four dimensional BF theory and gravity as shown in 

Two dimensional BF theory has been used as the basic theory in an attempt to define a manifold independent 
model of QFT in [23|. It is also related to grav ity in two dimensions in two ways: on the one hand it is equivalent to 
the so-called Jackiw-Teitelboim model [2TI l22j| . on the other hand it is related to usual 2d gravity via constraints in 
a way similar to the one exploited in four dimensions (see next section). The first relationship has been used in the 
canonical quantization of the Jackiw-Teitelboim model in 23] . The second relationship has been explored in (23 | 

Three dimensional BF theory and the spin foam quantization presented above is intimately related to classical and 
quantum gravity in three dimensions (for a classic reference see 25]). The state sum as presented above matches 
the quantum amplitudes first proposed by Ponzano and Regge in the 60 's based on their discovery of the asymptotic 
expressions of the 6j symbols and is often referred to as the Ponzano-Regge model. Divergences in the above 
formal expression require regularization. Natural regularizations are available and that the model is well defined 
[27I - I29I . For a detailed study of the divergence structure of the model see (sol - is^ . The quantum deformed version of 
the above amplitudes lead to the so called Turaev-Viro model [s^ which is expected to correspond to the quantization 
of three dimensional Riemannian gravity in the presence of a non vanishing positive cosmological constant. For the 
definition of observables in the latter context as well as in the analog four dimensional analog see 34] . 

The topological character of BF theory can be preserved by the coupling of the theory with topological defects that 
play the role of point particles. In the spin foam literature this has been considered form the canonical perspective in 
[35I [36| and from the covariant perspective extensively by Freidel and Louapre [st'I . These theories have been shown 
by Freidel and Livine to be dual, in a suitable sense, to certain non-commutative fields theories in three dimensions 

[ailli. 

Concerning coupling BF theory with non topological matter see [i^ |4l| for the case of fermionic matter, and [l^ 
for gauge fields. A more radical perspective for the definition of matter in 3d gravity is taken in [43| . For three 
dimensional supersymmetric BF theory models see [44l . |45| 

Recursion relations for the 6j vertex amplitudes have been investigated in (46l . |47| . They provide a tool for studying 
dynamics in spin foams of 3d gravity and might be useful in higher dimensions [48|. 



A. The coherent states representation 



In this section we introduce the coherent state representation of the SU (2) and Spin{4) path integral of BF theory. 
This will be particularly important for the definition of the models defined by Freidel and Krasnov in Q that we will 
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address in Section HVl as well as in the semiclassical analysis of the new models reported in Section IVlII The relevance 
of such representation for spin foams was first emphasized by Livine and Speziale in j49| . 

1. Coherent states 

Coherent states associated to the representation theory of a compact group have been studied by Thiemann and 
collaborators [SOySli. i51-59j see also [60| . Their importance for the new spin foam models was put forward by Livine 
and Speziale inpgjwhere the emphasis is put on coherent states of intertwiners or the so-called quantum tetrahedron 
(see also [filj). Here we follow the presentation of [63 . 

In order to built coherent states for Spin(4) we start by introducing them in the case of SU{2). Starting from the 
representation space Jifj of dimension dj = 2j + 1 one can write the resolution of the identity in tems of the canonical 
orthonormal basis |j, m) as 

Ij =^|j,m)(j, m|, (21) 

rn 

where —j<m< j. There exists an over complete basis \j,g) G J^j labehed by ^ G SU{2) such that 

lj =d, / dg\j,g){3,g\, (22) 



SU(2) 

The states |j, 17) S ^ are SU{2) coherent states defined by the action of the group on maximum weight states \j,j) 
(themselves coherent), namely 

\j,g)^g\j,j)^J2\^,m)Dl^{g), (23) 



m 

follows from the orthonormality of unitary representation matrix elements, namely 



where Dl^^{g) are the matrix elements of the unitary representations in the \j,m) (Wigner matrices). Equation 



d, / dg\j,g){],g\,^A,Y,\3.^){3.^'\ I c?5 (5)^^„', (5) = E 1^' 0' "^1' (^4) 
•^su(2) Jsu(2) 

where in the last equality we have used the orthonormality of the matrix elements. The decomposition of the identity 
p2|) can be expressed as an integral on the two-sphere of directions = SU {2) /U {!) by noticing that D^^ -^g) and 

D-'^-{gh) differ only by a phase for any group element h from a suitable U{1) C SU{2). Thus one has 

Ij ^dj / dn|j»(j,n|, (25) 

where n £ S*^ is integrated with the invariant measure of the sphere. The states |j, n) form (an over-complete) basis 
in ,3'ifj. SU{2) coherent states have the usual semiclassical properties. Indeed if one considers the generators of 
su{2) one has 

{j,n\.P\j,n)=jn\ (26) 

where is the corresponding three dimensional unit vector for n S^. The fluctuations of are also minimal with 
AJ^ = h^j, where we have restored h for clarity. The fluctuations go to zero in the limit h and j ^ 00 while 
Hj is kept constant. This kind of limit will be used often as a notion of semiclassical limit in spin foams. The state 
\j, n) is a semiclassical state describing a vector in R'^ of length j and of direction n. It will convenient to introduce 
the following graphical notation for equation ([25| 

j 

= d 



^0 



j dn n (27) 
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Finally, an important property of SU{2) coherent states stemming from the fact that 



L?'; j) 1272)1212)' 

is that 



1 i\ 



|j, = (28) 

The above property will be of key importance in constructing effective discrete actions for spin foam models. In 
particular, it will play a central role in the study of the semiclassical limit of the EPRLand FK modesl studied in 
Sections lllli and IIVI In the following subsection we provide an example for SpiniA) BF theory. 



2. SpiniA) BF theory: amplitudes in the coherent state basis 

Here we study the coherent states representation of the path integral for Spin{A) BF theory. The construction 
presented here can be extended to more general cases. The present case is however of particular importance for the 
study of gravity models presented in Sections IIIIl and IIVI With the introductions of coherent states one achieved the 
most difficult part of the work. In order to express the Spin{A) BF amplitude in the coherent state representation one 
simply inserts a resolution of the identity in the form (|25|) on each and every wire connecting neighbouring vertices 
in the expression (IT51) for the BF amplitudes. The result is 




where we have explicitly written the n± e S'^ integration variables only on a single cable. One observes that there 
is one n± G S'^ per each wire coming out at an edge e € A*; as wires are in one-to-one correspondence with faces 
/ G A* the integration variables rt^y G are labelled by an edge and face subindex. In order to get an expression of 
the BF path integral in terms of an affective action we restore at this stage the explicit group integrations represented 
by the boxes in the previous equation. One gets, 

zbaa)^ y: n d.jd^; / n d,.-d^-/-e7^n+ 

where we have used the coherent states property (|28p . and \n^) is a simplified notation for |i,n^). The previous 
equation can be finally written as 

Zbf{A)= Y 11^.7^,;/ n d,7/,+/"e7^</rf5e7rfffe'/ exp(5^,„±[5^]), (31) 
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where the discrete action 

^^.n±[5^]= E ^]:,n„[5^] (32) 
■uGA* 

with 

5 

^Inig]^ E 2ja6ln(nab|g^igb|nfc„), (33) 

a<b=l 

and the indices a, 6 label the five edges of a given vertex. The previous expression is exactly equal to the form (jlip 
of the BF amplitude. In the case of the gravity models studied in what follows, the coherent state path integral 
representation will be the basic tool for the study of the semiclassical limit of the models and the relationship with 
Regge discrete formulation of general relativity. 



B. The relationship between gravity and BF theory 

The field theory described in the present section has no local degrees of freedom. It represents the simplest example 
of a topological field theory in four dimensions. The interest of this theory for gravity model stems from the fact 
that an action for the gravitational degrees of freedom (basically equivalent to general relativity in the first order 
formulation) can be obtained by supplementing a 4d BF theory action with internal gauge group SL{2, C) (Lorentzian) 
or Spin{4:) (Riemannian) with the following set of quadratic constraints on the _B-field 

eijKLBliBfJ' - e e^,p. « 0, (34) 

where e = a^{l/Al)eijKLBj/l,BjfJ"e'^^P'^ where cr^ — ±1 according to we are in the Riemannian of Lorentzian case. 
More generally, a one parameter family of gravity actions can be obtained obtained from the imposition of the previous 
constraints on the following modified BF action 

S^{B,uj)^ f {CB+-B)AF{uj)), (35) 
Jm 7 

where 7 is the Immirzi parameter. The strategy behind the definition of the new spin foam models for quantum 
gravity consists of imposing these constraints on the path integral of BF theory on the momenta J = *B + ^B 
conjugate to w. In order to impose the Plebanski constraints above it will be convenient to express the B field in 
terms of the momenta J, namely 

i?=-^^(J-7V). (36) 

The imposition of the constraints (|67p on the BF path integral on a fixed discretization can be done in two different 
ways: by directly restricting the spin foam configurations (this is the EPRL approach described in the following 
section), or by restricting the semiclassical values of the B field in the coherent state representation of the BF path 
integral (this is the FK strategy described in Section HV)) . 



III. THE ENGLE-PEREIRA-ROVELLI-LIVINE (EPRL) MODEL 

In this section we introduce the Engle-Pereira-Rovelli-Livine (E PRL) model [3, The section is organized as 
follows. The relevant representation theory is introduced in Section UlI Al In Section UlIBI we present and discuss the 
linear simplicity constraints — classically equivalent to the Plebanski constraints — and discuss their implementation 
in the quantum theory. In Section IIII CI we introduce the EPRL model of Riemannian gravity. In Section IIII Dl we 
prove the validity of the quadratic Plebanski constraints — reducing BF theory to general relativity — directly in the 
spin foam representation. In Section IIII El we present the coherent state representation of the Riemannian EPRL 
model. In Section IIIIFI we describe the Lorentzian model. The material of this section will also allow us to describe 
the construction of the closely related (although derived from a different logic) Riemannian FK constructed in 
The idea that linear simplicity constraints are more convenient for dealing with the constraints that reduce BF theory 
to gravity was pointed out by Freidel and Krasnov in this last reference. 
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A. Representation theory of Spin{4) and SL{2, C) and the canonical basis 

In this section we present the representation theory of the groups Spin{4) and SL{2, C) that is neccesarry for the 
definition of the new spin foam models for Riemannian and Lorentzian gravity respectively. To emphasize the highly 
symmetric structure of the two we present them in a unified notation where a parameter a = 1 for the Riemannian 
sector and ct = i for the Lorentzian one. The simple relationship between the two might be a hint of a possible 
relationship between model amplitudes in a spirit similar to the interesting link between Euclidean and Lorentzian 
QFT provided by Wick rotations ^. Unitary irreducible representations J^p^k of Spin{A) and SL{2, C) are labelled by 
two parameters p and k. In the case of Spin{A) — SU{2) x SU{2) the unitary irreducible representations are finite 
dimensional and the labels p and k can be expressed in terms of the half integers labelling the right and left SU (2) 
unitary representations as follows 

P = ]++j- + l k=\j+-j-\. (37) 

In the SL{2, C) case the unitary irreducible representations arc infinite dimensional and one has 

peR+ fee N/2. (38) 

The two Casimirs are Ci = jJijJ^'^ = + a^K^ and C2 = j*JijJ^'^ = K ■ L where U are the generators of an 
arbitrary rotation subgroup and are the generators of the corresponding boosts. The Casimirs act on |p, k) G ^^p^k 
as follows 

Cib,fc)-i(fc2 + aV-l)b,A;) 

C2\p,k)^pk\p,k). (39) 



For detail on the representation theory of SL{2,C) see j63l - l65l |. The definition of the EPRL model requires the 
introduction of an (arbitrary) subgroup SU{2) C Spin{A) or SU{2) C SL{2, C) according to whether one is working in 
the Riemannian or in the Lorentzian sector. This subgroup corresponds to the internal gauge group of the gravitational 
phase space in connection variables in the time gauge (see Q for details). In the quantum theory, the representation 
theory of this SU(2) subgroup will be hence important. This importance will soon emerge as apparent from the 
imposition of the constraints that define the EPRL model. The link between the unitary representations of 5^(2, C) 
and those of SU{2) is given by the decomposition 

p-l j*+3^ 

'Km^®'^^3= (40) 

for the Riemannian sector, and 

00 

=^P,fc=0^-, (41) 
j=fc 

for the Lorentzian sector. As the unitary irreducible representations of the subgroup SU{2) £ Spin{4) and SU{2) G 
SL{2,C) are essential in understanding the link of the EPRL model and the operator canonical formulation of LQG 
it will be convenient to express the action of the generators of the Lie algebra of the corresponding group in a basis 
adapted to the above equation. In order to do this we first notice that the Lie algebra spin{4) and s?(2,C) can be 
characterized in terms of the generators of a rotation subgroup U and the remaining boost generators as follows 

[L3,L±] = ±L± [L+,L_]=2L3 
[L+,i^+] = [L-,i^-] = [L3,i^3] -0 

[ifa, i±] = ±K± [L±,K^] = ±2 i^a [^3, i^±] ^ ± K± 

[i^a, K±] = ±a^L± [K+,K^] = 2a' L^, (42) 



^ Such explicit relationship between gravity amplitudes in the Euclidean and Lorentzian sectors can be established by analytic continuation 
in 3d 1621. 
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where K± = ± iK"^ and L± = L^ztiL'^ respectively. The action of the previous generators in the basis \p, k;j, m) 
can be shown to be 

L^\p,k;j,m) = ■m\p,k;j,m), 

L+\p,k;j,m) = \/ {j + m + l)(j - m)\p, k; j, m + 1), 
L~\p,k;j,m} = ^ {j + m){j - m + l)\p, k; j,m ~ 1), 

K^\p, k;j,m) = a^^/f - m?\p, k;j - l,m) + j-jmlp, k;j,m) - aj+i\/(j + l)^ -JV?\p,k;j + l,m), 

K+\p,k;j,m) = [j - m){j - m - l)\p, k; j -l,m+l) 

+1] -m)(j + m + l)\p, k; j, m + 1) 

+aj+i\/ {j + m + + m + 2)\p, fc; j + 1, m + 1), 

K^\p,k;j,m) = -a^^ {j + m){j + m - l)\p,k;j - l,?7i- 1) 

+1] V + "i) {j -m+l)\p, k; j, m - 1) 

-ttj+i V(i -m+\){] -m + 2)\p,k-] + \,m-\), (43) 

where 



The previous equations will be important in what follows: they will allow for the characterisation of the solutions of 
the quantum simplicity constraints in both the Riemannian and Lorentzian models in a direct manner. This concludes 
the review of the representation theory that is necessary for the definition of the EPRL model. 



B. The linear simplicity constraints 



As first shown in 
presented below (p7|) - 



6|, the quadratic Plebanski simplicity constraints — and more precisely in their dual version 
-are equivalent in the discrete setting to the linear constraint on each face of a given tetrahedron 



1 



-K) « 0, 



(45) 



where the label / makes reference to a face / G A*, and where (very importantly) the subgroup SU{2) C Spin{4) 
or SL{2,C) that is necessary for the definition of the above constraints is chosen arbitrarily at each tetrahedron, 
equivalent on each edge e G A*. Such choice of the rotation subgroup is the precise analog of the time gauge in the 
canonical analysis of general relativity. The EPRL model is defined by imposing the previous constraints as operator 
equations on the Hilbert spaces defined by the unitary irreducible representations of the internal gauge group that 
take part in the state-sum of BF theory. We will show in Section IlII DI that the models constructed on the requirement 
of a suitable imposition of the linear constraints (|45l) satisfy the usual quadratic Plebanski constraints — that reduce 
BF theory to general relativity — in the path integral formulation (up to quantum corrections which are neglected in 
the usual semiclassical limit). 

From the commutation relations (I42p of previous section we can easily compute the commutator of the previous 
tetrahedron constraints and conclude that in fact it does not close, namely 



J 1 _ 



(1 + 7)^/ ~ 



1 



1 



7 ij jk 
2 — ^ k^f- 



(46) 



The previous commutation relations imply that the constraint algebra is not closed and cannot therefore be imposed 
as operator equations on the states summed over in the BF partition function in general. There are two interesting 
exceptions to the previous statement: 

1. The first one is to take 7 = ±(t. This corresponds to the description of the model in terms of self-dual or 
anti-self-dual variables. Unfortunately, the construction of the new models is not well defined in this case for 
the Lorentzian theory and leads to a trivial result in the Riemannian sector: SU{2) BF theory. 
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2. The second possibihty is to work in the sector where = 0. This choice leads to the Barret-Crane model [66[ 
where the degrees of freedom of BF theory seem over constrained: boundary states satisfying the BC constraints 
are a very small subset of the allowed boundary states in LQG. This is believed to be problematic if gravity is 
to be recovered at low energies. 

The EPRL model is obtained by restricting the representations appearing in the expression of the BF partition 
function so that at each tetrahedron the linear constraints (|45l) the strongest possible way that is compatible with 
the uncertainties relations stemming from (j46l) . In addition one would add the requirement that the state-space of 
tetrahedra is compatible with the state-space of the analogous excitation in the canonical context of LQG so that 
arbitrary states in the kinematical state of LQG have non trivial amplitudes in the model. 

Due to the fact that the constraints Dj do not form a closed (first class) algebra in the generic case one needs to 
devise a weaker sense in which they are to be imposed. One possibility is to consider the Gupta-Bleuler criterion 
consisting in selecting a suitable class of states for which the matrix elements on vanish. One notices from p3| 
that if we chose the subspace Jfj C J^p^k one has 

{p,k,j,q\Dl\p,k,j,m) = (5,,„im(l - ^) 

{p,k,j,q\Df\p,k,j,m) = Sq±l^„^^/{j ±ni + T m){l - ^). 
The matrix elements of the linear constraints vanish in this subclass if one can chose 

There are two cases: 

1. Case 7 < 1: Following [g^], in this case one restricts the representations to 

Riemannian: p = j + 1, k = jj. Lorentzian: p ~ j{j + 1), k = j. (48) 

which amounts to choosing the maximum weight component j = p — 1 in the expansion (|4T|) . In the Riemannian 
case the above choice translates into = (1 ± /2 for the SU(2) right and lect representations. Notice that 
the solutions to the simplicity constraints in the Riemannian and Lorentzian sectors look very different for 7 < 1. 
Simple algebra shows that condition (|47)) is met. There are indeed other solutions [11] of the Gupta-Blculcr 
criterion in this case. 

2. Case 7 > 1: In this case |69| one restricts the representations to 

Riemannian: p ~ ^{j + 1), fc = j. Lorentzian: p — ^{j + 1), k = j. (49) 

which amounts to choosing the minimum weight component j = fc in the expansion ()41|) and . For the Rieman- 
nian case we can write the solutions in terms of = (7 ± 1)| -I- Notice that for 7 > 1 there is complete 
symmetry between the solutions of the Riemannian and Lorentzian sectors. In my opinion, this symmetry 
deserves further investigation as it might be an indication of a deeper connexion between the Riemannian and 
Lorentzian models (again such relationship is a fact in 3d gravity |62j|). 

Another criterion for weak imposition can be developed by studying the spectrum of the Master constraint Mf = 
Df ■ Df. Strong imposition of the constraints would amount to looking for the kernel of the master constraint Mf . 
However, generically the positive operator associated with the master constraint does not contain the zero eigenvalue 
in the spectrum due to the open nature of the constraint algebra ()46|) . 

It is convenient 70] to express the master constraint in a in a manifestly invariant way. In order to get a gauge 
invariant constraint one starts from the master constraint and uses the i'j- — classically to write it in terms of 
Casimirs, namely 

Mf = (1 + C7^-/^)C2 - 2Ci7, 

where Ci and C2 are the Casimirs given in equation p9p . The minimum eigenvalue condition is 

Riemannian: p — j, k — ^j. Lorentzian: p — ■^j, k — j. (50) 

The minimum eigenvalue is m„i„ = h^jj{j'^ — 1) for the Riemannian case and m^j = 7 for the Lorentzian case. 
The master constraint criterion works better in the Lorentzian case as pointed out in I7OI • More recently, it has been 
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shown that the constraint solution p = "fj and k = j also follows naturally from a spinor formulation of the simplicity 
constraints (7ll - [73| . The above criterion is used in the definition of the EPRL model. 

It is important to point out that the Riemannian case imposes strong restrictions on the allowed values of the 
Immirzi parameter if one wants the spin j € N/2 to be arbitrary (in order to have all possible boundary states allowed 
in LQG). In this case the only possibilities are 7 = N or 7 = 1. This restriction is not natural from the viewpoint of 
LQG. Its relevance if any remains misterious at this stage. 

Summarising, in the Lorentzian (Riemannian) EPRL model one restricts the SL(2,<C) (S'pm(4)) representations of 
BF theory to those satisfying 



p^jj k = j 

for j G N/2. From now on we denote the subset of admissible representation 

^ C Irrep(5L(2,C))(Irrep(5pm(4))) 



(51) 



(52) 



The admissible quantum states ^ are elements of the subspace J^j C J^^-yjj (i.e., minimum weight states) satisfy the 
constraints P5)) in the following emiclassical sense: 



(53) 



where the symbol i^sc (order semiclassical) denotes a quantity that vanishes in limit ?i — > 0, j — > cxd with hj =constant. 
In the Riemannian case the previous equation can be written as 



[(i-7)^;-(i+7)^i]* = ^.c, 

which in turn has a simple graphical representation in terms of spin-network grasping operators, namely 



(54) 



-(I + 7) 



+ (1-7) 




(55) 




The previous equation will be of great importance in the graphical calculus that will allow us to show that the linear 
constraint imposed here at the level of states imply the vanishing of the quadratic Plebanski constraints (j67p and 
their fluctuations, computed in the path integral sense, in the appropriate large spin semiclassical limit. 



C. Presentation of the Riemannian EPRL amplitude 



Here we complete the definition of the EPRL models by imposing the linear constraints on the BF amplitudes 
constructed in Section (jlTI . We will also show that the path-integral expectation value of the Plebanski constraints 
([57]). as well as their fluctuations, vanish in a suitable semiclassical sense. This shows that the EPRL model can be 
considered as a lattice definition of the a quantum gravity theory. 

We start with the Riemannian model for which a straightforward graphical notation is available. The first step is 
the translation of equation (j40p — for p and k satisfying the simplicity constraints — in terms of the graphical notation 
introduced in Section (|lT]). Concretely, for 7 < 1 one has = (1 ± 7)j/2 G JC^ becomes 



(1-7)^ 



(1 + 7)! 



e 



(1-7) 




a 



,(l + 7)§ 



(56) 



(1 - i)l 



■Xl + 7)* 



For 7 > 1 we have 
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(7-1)! 



(l + 7)§ 



1] 



(7-1) 




(7-l')i 



a 



(1+7)! 



(57) 



tl+7)! 



The implementation of the hnear constraints of Section (IIIIBp consist in restricting the representations p/ of 
Spin(4) appearing in the state sum amphtudes of BF theory as written in Equation (jl8p to the subclass p/ £ J^-/ C 
Irrep(5'p«n(4)), defined above, while projecting to the highest weight term in (|56p for 7 < 1. For 7 > 1 one must 
take the minimum weight term in ([57| . The action of this projection will be denoted : ^^1+^)^/2, 1(1-7)1^/2 — ^ ; 
graphically 




(58) 



Explicitly, one takes the expression of the BF partition function (jl3p and modifies it by replacing the projector 
PinviPi^ • • • > P4) with pi, • • • p4 G ^ by a new object 

P!prl (Jl , • • • , J4) = Ptnv {pi---PiWn®---® % )PLu {Pl ' ' ' Pi) (59) 

with ji, ■ • ■ ji G N/2 implementing the linear constraints described in the previous section. Graphically the modifica- 
tion of BF theory that produces the EPRL model corresponds to the replacement 



PinviPl ' - ' P4) = HH 




on the expression (jlSp , where we have dropped the representation labels from the figure for simplicity. We have done 
the operation ()58p on each an every of the four pairs of representations. The Spin{4) integrations represented by 
the two boxes at the top and bottom of the previous graphical expression restore the full Spin(A) invariance as the 
projection ([55]) breaks this latter symmetry for being based on the selection of a special subgroup SU{2) C Spin{4) 
in its definition (see section ?? for an important implication). One should simply keep in mind that green wires in 
the previous two equations and in what follows are labeled by arbitrary spins j (which are being summed over in 
the expression of the amplitude (l6T|) ). while red and blue wires are labelled by j+ = (1 + 7)j/2 and = |1 — 7|j72 
respectively. With this ([TS]) is modified to 

pjex feA* e 




The previous expression is defines the EPRL model amplitude. 
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1. The spin foam representation of the EPRL amplitude 

Now we will work out the spin foam representation of the EPRL amplitude which at this stage will take no much 
more effort than the derivation of the spin foam representation for Spin{'i) BF theory as we went from equation (|18p 
to (Uni) in Section im The first step is given in the following equation 




which follows basically from the invariance of the Haar measure ^ (in the last line we have used (ITTl) ) . More presicely, 
the integration of the subgroup SU{2) G Spin{4), represented by the green box on the right, can be absorbed by 
suitable redefinition of the integration on the right and left copies of SU{2), represented by the red and blue boxes 
respectively. With this we can already write the spin foam representation of the EPRL model, namely 




where the vertex amplitude (graphically represented) depends on the 10 spins j associated to the face-wires and the 
5 intertwiners associated to the five edges (tetrahedra). As in previous equations we have left the spin labels of wires 
implicit for notational simplicity. We can write the previous spin foam amplitude in another form by integrating out 
all the projectors (boxes) explicitly. Using, ([TTl) we get 
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thus replacing this in (pT|) we get 

jf feA* te v£A* 




where the coefficients /^V^- are the so-caUed fusion coefficients which appear in their graphical form already in (|64p . 
more explicitly 




The previous is the form of the EPRL model as derived in 

D. Proof of validity of the Plebanski constraints 

In this section we prove that the quadratic constraints are satisfied in the sense that their path integral expectation 
value and fluctuation vanish in the appropriate semiclassical limit. 

1. The quadratic Plebanski constraints 
The quadratic Plebanski constraints are 

<^iJKLBjJ,BpJ^ - e e^^pcr w 0. (67) 
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The constraints in this form are more suitable for the translation into the discrete formulation. More precisely, 
according to ©, the smooth fields Bj^'l is now associated with the discrete quantities B^J^^^^^^^, or equivalently By as, 
we recall, faces / S A"^ are in one-to-one correspondence to triangles in four dimensions. The constraints ([67]) are local 
constraints valid at every spacetime point. In the discrete setting, spacetime points are represented by four-simplexes 
or (more addapted to our discussion) vertices v £ A*. With all this the constraints ((67)) are discretized as follows: 

Triangle (or diagonal) constraints: ^ijklBY Bj^ = 0, (68) 

for all / e w, i.e., for each and every face of the 10 possible faces touching the vertex v. 

Tetrahedron constraints: ^ijKLBy B^,^ ~ 0, (69) 

for all /, /' G u such that they are dual to triangls sharing a one-simplex, i.e., belonging to the same tetrahedron out 
of the five possible ones. 

4-simplex constraints: ^ijklBY B^^ — e„, (70) 

for any pair of faces f,fGv that are dual to triangles sharing a single point. The last constraint will require a more 
detailed discussion. At this point let us point out that the constraint (j70p is interpreted as a definition of the four 
volume By of the four-simplex. The constraint requires that such definition be consistent, i.e., the true condition is 

^ijklBY B^^ — eijKhBy B^,^ — eijKLBj'lBp^ — ■ ■ ■ — (71) 

for all five different possible pairs of / and / in a four simplex, and where we assume the pairs /-/ are ordered in 
agreement with the orientation of the complex A*. 

2. The path integral expectation value of the Plebanski constraints 

Here we prove that the Plebanski constraint are satisfied by the EPRL amplitudes in the path integral expectation 
value sense. 

The tnangle constraints: 

We start from the simplest case: the triangle (or diagonal) constraints We choose a face f € v (dual 

to a triangle) in the cable- wire-diagram of Equation (j61l) . This amounts to choosing a pair of wires (right and 
left representations) connecting two nodes in the vertex cable wire diagram. The two nodes are dual to the two 
tetrahedra — in the four simplex dual to the vertex — sharing the chosen triangle. From equation (j36p can show that 

euKLB'/Bf"^ oc (14- 7)^7 • Jy - (1 - 7)^+ • J+, (72) 

where Jj^ denotes the self-dual and anti-self-dual parts of Tiy . The path integral expectation value of the triangle 
constraint is then 

((1 + 7)V7. J--(l-7)2j+- J+)cx (73) 




where the double graspings on the anti-self-dual (blue) wire and the self-dual (red) wire represent the action of the 
Casimirs J J ■ J J and • on the cable- wire diagram of the corresponding vertex. Direct evaluation shows that the 
previous diagram is proportional to fi^jf which vanishes in the semiclassical limit h 0, j ^ 00 with hj ^constant. 
We use the notation already adopted in ([5^ and call such quantity ^Tgc- This concludes the proof that the triangle 
Plebanski constraints are satisfied in the semiclassical sense. 
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The tetrahedra constraints: 
The proof of the vahdity of the tetrahedra constraints (|69)) . In this case we also have 




where we have chosen an arbitrary pair of faces. In order to prove this let us develop the term on the right. The 
result follows from 




where in the first line we have used the fact that the double grasping can be shifted through the group integration 
(due to gauge invariance ([9|)), and in the first and second terms on the second line we have used Equation ([55|) to 
move the graspings on self-dual wires to the corresponding anti-self-dual wires. Equation ()74|) follows immediately 
from the pervious one; the argument works in the same way for any other pair of faces. Notice that the first equality 
in Equation ((75| implies that we can view the Plebanski constraint as applied in the frame of the tetrahedron as 
well as in a Lorentz invariant framework (the double grasping defines an intertwiner operator commuting with the 
projection Pf^^ represented by the box). An analogous statement also holds for the triangle constraints ([75]) . 

The 4-stmplex constraints 

Now we show the validity of the four simplex constraints in their form (|7ip . As we show below, this last set of 
constraints follow from the Spin{4) gauge invariance of the EPRL node (i.e., the validity of the Gauss law) plus the 
validity of the tetrahedra constraints ((69|) . Gauge invariance of the node takes the following form in graphical notation 
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where the above equation represents the gauge invariance under infinitesimal left SU{2) rotations. An analogous 
equation with insertions on the right is also valid. The validity of the previous equation can again be related to the 
invariance of the Haar measure used in the integration on the gauge group that defines the boxes (j9]). 

Now we chose an arbitrary pair / and / (where, recall, / is one of the three possible faces whose dual triangle 
only shares a point with the corresponding to /) and will show how the four volumen defined by it equals the one 
defined by any other admissible pair. The first step is to show that wc get the same result using the pair /-/ and 
/-/, where / is another of the three admissible faces opposite to /. The full result follows from applying the same 
procedure iteratively to reach any admissible pair. It will be obvious from the treatment given below that this is 
possible. Thus, for a given pair of admissible faces we have 




where going from the first line to the second and third lines we have simply used (I76p on the bottom graspings on 
the right and left wires. The last line results from the validity of (|69|): notice that the second terms in the second 
and third lines add up to as well as the third terms in the second and third line. There is an overall minus sign 
which amounts for an orientation factor. It should be clear that we can apply the same procedure to arrive at any 
admissible pair. 
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3. Pepri is not a projector 

Let us study in a bit more detail the object P^p^i (ii , • • • , J4)- We see that it is made of two ingredients. The first one 
is the projection to the maximum weight subspace J^j for 7 > 1 in the decomposition of J^j+j- for = (1 ± j)j/2 
(j^ = (7±l)j'/2 for 7 > 1) in terms of irreducible representations of an arbitrarily chosen SU{2) subgroup of Spin{A). 
The second ingredient is to eliminate the dependence on the choice of subgroup by group averaging with respect to 
the full gauge group Spin{4:). This is diagramaticaly represented in (|60l) . However P^p^iiji, ■ ■ • , J4) is not a projector, 
namely 

P!prl (Jl , • • • , J4)' ^ P!prl (jl , • • • , H)- (78) 

Technically this follows from ([5^ and the fact that 

[^;.(Pi---P4),(%®---®^J]^0 (79) 

i.e., the projection imposing the linear constraints (defined on the frame of a tetrahedrom or edge) and the S'pm(4) 
(or Lorentz) group averaging — rendering the result gauge invariant — do not commute. The fact the P^p^iiji, ■ ■ • ,34) 
is not a projection operator has important consequences in the mathematical structure of the model: 

1. From ([6T|) one can immediately obtain the following expression for the EPRL amplitude 

Zepw(A)= ^ n d|l-^||Cl(i+^)|n^4w(jl'--- '^4). (80) 

This expression has the formal structure of expression (|13l) for BF theory. The formal similarity however is 
broken by the fact that PepriUi^ ' ' ' 7J4) is not a projection operator. From the formal perspective is the 
possibility that the amplitudes be defined in term of a network of projectors (as in BF theory) might provide 
an interesting structure that might be of relevance in the definition of a discretization independent model. On 
the contrary, the failure of P^priUij ' " ' iji) to be a projector may lead, in my opinion, to difficulties in the limit 
where the complex A is refined: the increasing of the number of edges might produce either trivial or divergent 
amplitudes ^. 

2. Another difficulty associated with P^p^iiji, ■ ■ ■ , J4)^ 7^ Pepriij^^ ' ' ' 'i*) failure of the amplitudes of the 
EPRL model, as defined here, to be consistent with the abstract notion of spin foams as defined in [t^I- This is a 
point of crucial importance under current discussion in the community. The point is that the cellular decomposi- 
tion A has no physical meaning and is to be interpreted as a subsidiary regulating structure to be removed when 
computing physical quantities. Spin foams configurationa can fit in different ways on a given A, yet any of these 
different embeddings represent the same physical process (like the same gravitational field in different coordi- 
nates). Consistency requires the spin foam amplitudes to be independent of the embedding, i.e., well defined on 
the equivalence classes of spin foams as defined by Baez in (t^ I (the importance of these consistency requirements 
was emphasized in [ill)- The amplitude (15(11) fails this requirement due to PLriiJii " ' 7^ Pepriih: ' ■ • iJi)- 



4- The Warsaw proposal 

If one sees the above as difficulties then there is a simple solution, at least in the Riemannian case. As proposed in 
[tI,!?^! one can obtain a consistent modification of the EPRL model by replacing P^p^i in (|5D|) by a genuine projector 
Pw, graphically 




^ This is obviously not clear from the form of 1 I8OII . We are extrapolating the properties of (Pf ^i)^ ^'-'^ large to those of the amplitude 
II8OI I in the large number of edges limit implied by the continuum limit. 
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It is easy to check that by construction 

The variant of the EPRL model proposed in [tI, [t^I takes then the form 

Zepw(A) d|i_,|.d(i_,^).n^^(j'l'---'^4) 



(82) 



(83) 



EE n 



-711^(1+7)^ 



n n 



3f 



/GA* 



eeA* 



uGA* 




Thus in the modified EPRL model edges e G A* are assigned pairs of intertwiner quantum numbers and and 
an edge amplitude given by the matrix elements gf^ (where Vs and Vt stand for the source and target vertices of 

the given oriented edge). The fact that edges are not assigned a single quantum number is not really significative; 
one could go to a basis of normalized eigenstates of and rewrite the modified model above as a spin foam model 
where edges are assigned a single (basis element) quantum number. As the nature of such basis and the quantum 
geometric interpretation of its elements is not clear at this stage, it seems simpler to represent the amplitudes of the 
modified model in the above form. 

The advantages of the modified model are important,; however, a generalization of the above modification of the 
EPRL model in the Lorentzian case is still lacking. Notice that this modification does not interfere with the results on 
the semiclassical limit (to leading order) as reviewed in Section IVIII The reason is that the matrix elements Sa/s 
in that limit ItsI. 



E. The coherent states representation 



We have written the amplitude defining the EPRL model by constraining the state sum of BF theory. For semi- 
classical studies that we will review in Section IVIII it is convenient to express the EPRL amplitude in terms of the 
coherent states basis. The importance of coherent states in spin foam models was put forward in [49!] and explicitly 
used to re-derive the EPRL model in [t^. The coherent state technology was used by Freidel and Krasnov in to 
introduce a new kind of spin foam models for gravity: the FK models. In some cases the FK model is equivalent to 
the EPRL model; we will review this in detail in Section lIVl 

The coherent state representation of the EPRL model is obtained by replacing (|27p in each of the intermediate 
SU{2) (green) wires in the expression (|5T|) of the EPRL amplitudes, namely 




where we used in the last hne the fact that for 7 < 1 the representations j of the subgroup SU{2) G Spin{A) are 
maximum weight, i.e., j = + . Doing this at each edge we get 
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3f feA* 




where we have exphcitly written the n G integration variables on a single cable. The expression above is very 
similar to the coherent states representation of Spin(4:) BF theory given in Equation (|29p . In fact one would get 
the above expression if one would start form the expression (j29p and would set nj^ = = n^f while dropping 

for example all the sphere integrations corresponding to the (or equivalently n'^f)- Moreover, by construction 
the coherent states participating in the previous amplitude satisfy the linear constraints (pS)) in expectation values, 
namely 

{j,nef\D}\j,nef) = 

= (j, - 7) + (1 + 7) J/Ij, n,f) = 0. (87) 

Thus the coherent states participating in the above representation of the EPRL amplitudes solve the linear simplicity 
constraints in the usual semiclassical sense. The same manipulations leading to (f89|) in Section HI] lead to a discrete 
effective action for the EPRL model, namely 

Kri-T. n d(,_^)i,d^,^^^i, / n ^J.fdnefdg:fdgtj exp {S]^^Jg^]) , (88) 

3f feA* ^ ^ eeA* 

where the discrete action 

vGA* 

with 

5 

Sln[9]= E 2jab In (nab I ffa^ 56 1 «6a), (90) 
a<b=l 

and the indices a, b label the five edges of a given vertex. The previous expression is exactly equal to the form (fTTI) of 
the BF amplitude. In the case of the gravity models presented here, the coherent state path integral representation 
(analogous to ((3T|) ) will be the basic tool for the study of the semiclassical limit of the models and the relationship 
with Regge discrete formulation of general relativity. 

The case 7 > 1 

The case 7 > 1 is more complicated [s^. The reason is that the step (|85]) directly leading to the discrete action 
in the previous case is no longer valid as the representations of the subgroup SU(2) G Spin{A) are now minimum 
instead of maximum weight. However, the representations j+ = + j are maximum weight. We can therefore insert 
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where we are representing the relevant part of the diagram appearing in equation (I84|) . In the last line we have used 
that j+ = j + (i.e. maximum weight), and the graphical notation mo-o n = (^rn\n) as it follows from our previous 
conventions. With all this, one gets 

Zi;r\ = (92) 
E n d(,_^)£xd(i+^)^ / n ^J.f\i+^y^'in-fdmefdg;fdg+ exp (^;±>;„ „^[g±]), 

3f /GA* eeA* 

where the discrete action 

S]^:r.j9^] = E S]^,n,m[9^] (93) 

with 

S]±,n,m{9^] = E jab{l + l)l0g{{mab\gab\'mba)) + jabil ~ 1) log( (Wafc Iff^b ) ) + 

l<a<h<5 

+ 2jab {log{{Tlab\mab)) + log( (to6q l^ba) ) ) • (94) 



1. Some additional remarks 



It is important to point out that the commutation relations of basic fields — reflecting the simple algebraic structure 
of spin{4) — used here is the one induced by the canonical analysis of BF theory presented previously. The presence of 
constraints generally modifies canonical commutation relations in particular in the presence of second class constraints. 
For some investigation of the issue in the context of the EPRL and FK models see [69i| . In "sT] it is pointed out that 
the presence of secondary constraints in the canonical analysis of Plebanski action should translated in additional 
constraints in the holonomies of the spin foam models here considered (see also (s^l). A possible view is that the 
simplicity constraints are here imposed for all times and thus secondary constraints should be imposed automatically. 

There are alternative derivations of the models presented in the previous sections. In particular one can derive 
them from a strict Lagrangean approach of Plebanski's action. Such viewpoint is taken in [83l - l85| . The path integral 
formulation of Plebansky theory using commuting _B-fields was studied in [sB], where it is shown that only in the 
appropriate semiclassical limit the amplitudes coincide with the ones presented in the previous sections (this is just 
another indication that the construction of the models have a certain semiclassical input; see below). The spin foam 
quantization of the Hoist formulation of gravity via cubulations was investigated in [87| . The simplicity constraints 
can also be studied from the perspective of the U{N) formulation of quantum geometry Such U{N) treatment 
is related to previous work [89j, |90| which has been extended to a completely new perspective on quantum geometry 
with possible advantageous features [9l|, [o^] . For additional discussion on the simplicity constraints see [93,] . 
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Presentation of the EPRL Lorentzian model 



As briefly discussed in Section IIII[ unitary irreducible representations of 5i(2,C) are infinite dimensional and 
labelled by a positive real number p G and a half-integer k G N/2. These representation are the ones that 
intervene in the harmonic analysis of square integrable functions of SL(2,C) [6^ . Consequently, one has an explicit 
expression of the delta function distribution (defined on such test function) , namely 



jmjni 



(95) 



where Dj^j,^,{g) with j > k and j > m > — j (similarly for the primed indices) are the matrix elements of the 
unitary representations p — k in the so-called canonical basis [63| . One can use the previous expression the Lorentzian 
version of Equation (jlip in order to introduce a formal definition of the BF amplitudes, which now would involve 
integration of the continuous labels p/ in addition of sums over discrete quantum numbers such as k, j and ni. The 
Lorentzian version of the EPRL model can be obtained from the imposition of the linear simplicity constraints to 
this formal expression. As the continuum labels pf are restricted to pf = "fjf the Lorentzian EPRL model becomes 
a state-sum model as its Riemannian relative. Using the following graphical notation 



D 



jnij'm 



,{g)= l,rn 
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(96) 



the amplitude is 



3f /6A 




where the boxes now represent SL{2, C) integrations with the invariant measure. The previous amplitude is equivalent 
to the its spin foam representation 



^fpw(A)=EE n (1+7^)^1 n 



/SA* 



■ugA* 




The vertex amplitude is divergent due to the presence of a redundant integration over 5'£(2, C), it becomes finite by 
dropping an arbitrary integration, i.e. removing any of the 5 boxes in the vertex expression [94| . 



1. The coherent state representation 



It is immediate to obtain the coherent states representation of the Lorentzian models. As in the Riemannian case, 
one simply inserts resolution of the identities ([^^ on the intermediate SU{2) (green) wires in Wf\ from where it 
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results 



Jf /6A* 



n 




(97) 



IV. THE FREIDEL-KRASNOV (FK) MODEL 



Shortly after the appearance of the paper [4j , Freidel and Krasnov ^6)] introduced a set of new spin foam models for 
four dimensional gravity using the coherent state basis of the quantum tetrahedron of Livine and Speziale [1^ . The 
idea is to impose the linearized simplicity constraints (|45p directly as a semiclassical condition on the coherent state 
basis. As we have seen above, coherent states are quantum states of the right and left tetrahedra in BF theory which 
have a clear-cut semiclassical interpretation through their property We have also seen that the imposition of 

the linear constraints (|45l) a la EPRL is in essence semiclassical as they are strictly valid only in the large spin limit. 
In the FK approach one simply accept from the starting point that, due to their property of not defining set that 
is closed under commutation relations, the Plebansky are to be imposed semiclassically. One defines new models by 
restricting the set of coherent states entering in the coherent state representation of Spin{4) BF theory p9|) to those 
that satisfy condition (pSll in expectation values. They also emphasize how the model [1] corresponds indeed to the 
sector 'f — oo which has been shown to be topological [95| . 



The case 7 < 1 



For 7 < 1 the vertex amplitude is identical to the EPRL model. This is apparent in the coherent state expression 
of the EPRL model (EHl). Thus we have 
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(98) 



From the previous expression we conclude that the vertex amplitudes of the FK and EPRL model coincide for 7 < 1 

(99) 
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Notice however that different weights are assigned to edges in the FK model. This is due to the fact that one is 
restricting the Spin{4:) resolution of identity in the coherent basis in the previous expression, while in the EPRL 
model the coherent state resolution of the identity is used for SU{2) representations. This difference is important and 
has to do with the still un-settled discussion concerning the measure in the path integral representation. 

The case 7 > 1 

For the case 7 > 1 the FK amplitude is given by 

jf feA* 
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eSA* 



(l+7)i^(^-l)^ 
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m 



(100) 



The study of the coherent state representation of the FK model for 7 > 1, and comparison with equation (PT|) for the 
EPRL model, clearly shows the difference between the two models in this regime. 

Zh-H n d(,_^)^d^^^^^^ / n 'i\,.,\l^\,^,y:^dn^fd9:fdgtf exp (5/1 (101) 

jf /GA* •' eGA* 

where the discrete action 

^■£;[5^]= E(^n lal + S^^^., ^ \g+]), (102) 

J >" ^ (1-7)-^ ,n (l+7)^,s(7)n 

where 5(7) = sign(l — 7) and 

5 

S]^n[9]^ E 2jaf,ln(naf,|ffQ ^5f)|"-f)a), (103) 
a<b=l 

with the indices a, b labelling the five edges of a given vertex. 



V. BOUNDARY DATA FOR THE NEW MODELS AND RELATIONSHIP WITH THE CANONICAL 

THEORY 

So far we have considered cellular complexes with no boundary. Transition amplitudes are expected to be related to 
the definition of the physical scalar product. In order to define them one needs to consider complexes with boundaries. 
Boundary states are defined on the boundary of the dual two-complex A* that we denote dA*. The object dA* is a 
one-complex (a graph). According to the construction of the model (Section IIIII) boundary states are in one-to-one 
correspondence with SU (2) spin networks. This comes simply from the fact that links (one-cells) i G dA* inherit the 
spins labels (unitary irreducible representations of the subgroup SU{2)) of the boundary faces while nodes (zero-cells) 
n G dA* inherit the intertwiner levels of boundary edges. 

At this stage one can associate the boundary data with elements of a Hilbert space. Being in one-to-one correspon- 
dence with SU (2) spin networks, a natural possibility is to associate to them an element of the kinematical Hilbert 
space of LQG. More precisely, with a given coloured boundary graph 7 with links labelled by spins ji and nodes 
labelled by interwiners /,„ we associate a cylindrical function ^j.{jt},{i.„} G [SU [2)^ ^) , where here iV^ denotes 
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number of links in the graph 7. In this way, the boundary Hilbert space associated with 9 A* is isomorphic (if one 
used the natural AL measure) with the Hilbert space of LQG truncated to that fixed graph. Moreover, geometric 
operators such as volume and area defined in the covariant context are shown to coincide with the corresponding 
operators defined in the canonical formulation [6^, • Now, if cellular complexes are dual to triangulations then 
the boundary spin networks can have at most four valent nodes. This limitation can be easily overcome: as in BF 
theory the EPRL amplitudes can be generalized to arbitrary complexes with boundaries given by graphs with nodes 
of arbitrary valence. The extension of the model to arbitrary complexes has been first studied in |97l . |98| , it has also 
been revisited in [68j . 

Alternatively, one can associate the boundary states with elements oi ^^{Spin{4)^^) (in the Riemannian models) — 
or carefully define the analog of spin network states as distributions in the Lorentzian case (see [99*1 for some insights 
on the problem of defining a gauge invariant Hilbert space of graphs for non compact gauge groups). In this case 
one gets special kind of spin network states that are a subclass of the so-called projected spin networks introduced in 
Il0l| in order to define an heuristic quantization of the (non -con i muta.tive and very complicated) Dirac algebra 
of a Lorentz connection formulation of the phase space of gravity [lOCt Il02l - ll07l | . The fact that the se sp ecial subclass 
of projected spin networks appear naturally as boundary states of the new spin foams is shown in [108|. 

Due to their similarity for 7 < 1 the same relationship between boundary data and elements of the kinematical 
Hilbert space hold for the FK model. However, the such simple relationship does not hold for the model in the case 
7 > 1. 

It is important to mention that the knot ting properties of boundary spin network do not seem to play a role in 
present definitions of transition amplitudes |1Q9| | . 



VI. FURTHER DEVELOPMENTS AND RELATED MODELS 



The spin foam amplitudes discussed in the previous sections have been introduced by constraining the BF histories 
through the simplicity constraints. However, in the path integral formulation, the presence of constraints has the 
additional effect of modifying the weights with which those histories are to be summed: second class constraints 
modify the path integral measure (in the spin foam context this issue was raised in [zHi)- As pointed out before, this 
question has not been completely settled yet in the spin foam community. The explicit modification of the formal 
measure in terms of continuous variables for the Plebans ky ac tion was p resented in IIQJ. A systematic investigation 
of the measure in the spin foam context was attempted in and |112| . As pointed out in [T^ : there are restrictions 
in the manifold of possibilities coming from the requirement of background independence. The simple BF measure 
chosen in the presentation of the amplitudes in the previous sections satisfy these requirements. There are other 
consistent possibilities; see for instance |ll3l | for a modified measure which remains extremely simple and is suggested 
from the structure of LQG. 

An important question is the relationship between the spin foam amplitudes and the canonical operator formulation. 
The question of wether one can reconstruct the Hamiltonian constraints out of spin foam amplitudes has been ana lysed 
in detail in three dimensions. For the study of quantum three dimensional gravity from the BF perspective see [HJ], 
we will in fact present this perspective in detail in the three dimensional part of this article. For the relat ionship 
with the canonical theory using variables that are natural from the Regge gravity perspective see [llSl Ill6| There 
are generalizations of Regge variables moro adapted to the interpretation of spin foams [11 7J . In four dimensions the 
question has been investi gated in [118| in the context of the new spin foam models. In the context of group field theories 
this issue is explored in [ll9|. Finally, spin foams can in principle be obtained directly from the implementation of 
the Dirac program using path inte gral methods this avenue has been explored in [120, 121] from which a discrete path 
integral formulation followed [l22j | . The questi on o f the relationship between covariant and canonical formulations in 
the discrete setting has been analyzed also in [12^ . 

By construction all tetrahedra in the FK and EPRL models are embedded in a spacelike hypersurface and hence 
have only spacelike triangles. It seem natural to ask the question of whether a more general construction allowing 
for timelike faces is possible. The models described in previous sections have been generalized in order to include 
timelike faces in the work of F. Conr adv jl24| - [l26| . An earlier attempt to define such models in the context of the 
Barrett-Crane model can be found in |l27l |. 

The issue of the coupling of the new spin foam models to matter remains to a large extend un-explored territory. 
Nevertheless some results can be found in the literature. The c oupling of the Barrett-Crane model (the 7 — > cxd limit 
of the EPRL model) to Yang-Mills fields was studied in (l28j |. More recently the coupling of the EPRL model to 
fermions has been investi gate d in [129, 130]. A novel possibility of unification of the gravitational and gauge fields 
was recently proposed in jl3lj . 

The introduction of a cosmological constant in the construction of four dimensional spin foam models has a long 
history. Barrett and Crane introduced a vertex amplitude jl32| ] in terms of the Crane and Yetter model [l^ for BF 
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theo ry with cosmological constant. The Lorentzian quantum deformed version of the previous model was studied in 
[l33j . For the new models the coupli ng w i th a cosmological constant is explore d in terms of the quantum deformation 
of the internal gauge symmetry in |l34l . Il35l | as well as (independently) in |l36l |. The asymptotic s of the vertex 
amplitude are shown to be consistent with a cosmological constant term in the semiclassical limit in [l37'|. 

The spin foam approach applied to quantum cosmology has been explored in ^138r,143j . Th e spin foam formulation 
can also be obtained from the canonical pi cture prov ided by loop quantum cosmology (see |l44l | and references therein) . 
This has been explored systematically in [I45l4l48| . 

As we have discussed in the introduction of the new models, Heisenberg uncertainty principle precludes the strong 
imposition of the Plebanski constraints that reduce BF theory to general relativity. The results on the semiclassical 
limit of these models seem to indicate that metric gravity should be recovered in the low energy limit. However, its 
seerns lik ely that the semiclassical limit could be related to certain modifications of Plebanski's formulation of gr avity 
' l49l4l53 [. A simple interpretation of the new models in the context of the bi-gravity paradigm proposed in |l54l | 



could be of interest. 

As already pointed out in [t^ ] spin foams can be interpreted in close analogy to Feynman diagrams. Standard 
Fcynman graphs are generalized to 2-complexes and the labeling of propagators by momenta to the assignment 
of spins to faces. Finally, momentum conservation at vertices in standard feynmanology is now represented by 
spin-conservation at edges, ensured by the assignment of the corresponding intertwiners. In spin foam models the 
non-trivial content of amplitudes is contained in the vertex amplitude which in the language of Feynman diagrams 
can be interpreted as an intera ction. Th is analogy is indeed realized in the formulation of spin foam models in terms 
of a group field theory (GFT) [T55l . [T56j . 

The GFT formulation resolves by definition the two fundamental conceptual problems of the spin foam approach: 
diffeomorphism gauge symmetry and discretization dependence. The difficulties are shifted to the question of the 
physical role of A and the convergence of the corresponding perturbative series. 

In three dimensions this idea has been studied in more detail. In (l57j | scaling properties of the modification of the 
Boulatov group field theory introdu ced in |l58l | was studied in detail. In a further modification of the previous model 
(known as coloured tensor models [159| ') new techniques based on a suitable 1/7V expansion imply that amplitudes 
are dominated by spheric al topolo gy [160j; moreover, it seem possible that the continuum limit might be critical as 
in certain matrix models jl6ll4l65 |. However, it is not yet clear if there is a sense in which these models correspond 
to a physical theory. The naive interpretation of the models is that they correspond to a formulation of 3d quantum 
gravity including a dynamical topology. 



VII. RESULTS ON THE SEMICLASSICAL LIMIT OF EPRL-FK MODELS 



Having introduced the relevant spin foam models in the previous sections we now present the results on the large 
spin asymptotics of the spin foam amplitudes suggesting that on a fixed discretiza tion the semiclassical limit of the 
EPRL-FK models is given by Regge's discrete formulation of general relativity [sol Il66j |. 

The semiclassical limit of spin foams is based on the study of the the large spin limit asymptotic behaviour of 
coherent state spin foam amplitudes. The notion of large spin can be defined by the rescaling of quantum numbers 
and Planck constant according to j — > Xj and h — > h/X and taking A >> 1. In this limit the quantum geometry 
approximates the classical one when tested with suitable states (e.g. coherent states). However, the geometry remains 
discrete during this limiting process as the limit is taken on a fixed regulating cellular structure. That is why one 
usually makes a clear distinction between semiclassical limit and the continuum limit. In the semiclassical analysis 
presented here one can only hope to make contact with discrete formulations of classical gravity; hence the importance 
of Regge calculus in the discussion of this section. 

The key technical ingredient in this analysis is the repre sentation of spin foam amplitudes in terms of the coherent 
state basis introduced in Section III Al Here we follow |80l. [l66l - [T69l |. The idea of using coherent states and discrete 
effective actions for the study of the large spin asymptotics of spin foam amplitudes was put forward in (l70l Il7lj . The 
study of the large spin asymptotics has a long tradition in the context of quantum gravi ty dating back to the studied 
of Ponzano- Regge [2y| . More directly related to our discussion here are the early works [l72l Il73j | . The key idea is to 
use asymptotic stationary phase methods for the amplitudes written in terms of the discrete actions presented in the 
previous section. 

In this section we review the results of the analysis of the large spin asymptotics of the EPRL vertex amplitude for 
both the Riemannian and Lorentztian models. We follow the notation and terminology of 80] and related papers. 
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2. SU(2) 15j-symbol asymptotic 



As SU{2) BF theory is quite relevant for the construction of the EPRL-FK models, the study of the large spin 
asymptotics of the SU{2) vertex amplitude is a key ingredient in the analysis of [80]. The coherent state vertex 
amplitude is 



15j{j,n)= Yl^9a n (^ab\ga^9b\nbaf''' 



(104) 



a=l l<a<&<5 

which depends on 10 spins jab and 20 normals Uab ^ nba- The previous amplitude can be expressed as 
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and the indices a, h label the five edges of a given vertex. The previous expression is exactly equal to the form (jlip 
of the BF amplitude. In the case of the EPRL model studied in Sections IIIII the coherent state representation — see 
equations |88l |92l and [97] — is the basic tool for the study of the semiclassical limit of the models and the relationship 
with Regge discrete formulation of general relativity. 

In order to study the asymptot ics of (11051) one needs to use extended stationary phase methods due to the fact the 
the action (|106l) is complex (see |l70l . Il7l| ). The basic idea is that in addition to stationarity one requires real part 
of the action to be maximal. Points satisfying these two conditions are called critical points. As the real part of the 
action is negative definite, the action at critical points is purely imaginary. 

Notice that the action (I106P depends parametrically on 10 spins j and 20 normals n. These parameters define the 
so-called boundary data for the four simplex v G A*. Thus, there is an action principle for every given boundary 
data. The number of critical points and their properties depend on these boundary data, hence the asymptotics of 
the vertex amplitude is a function of the boundary data. Different cases are studied in detail in (80| . here we present 
their results in the special case where the boundary data describe a non-degenerate Regge geometry for the boundary 
of a four simplex, these data are referred to as Regge-like, and satisfy the gluing constraints. For such boundary data 
the action ()106|) has exactly two critical points leading to the asymptotic formula 
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(107) 



where Qab the appropriate diahedral angles defined by the four simplex geometry; finally the N± are constants that 
do not scale with A. 

3. The Riemannian EPRL vertex asymptotics 

The previous result together with the fact that the EPRL amplitude for 7 < 1 is a product of SU (2) amplitudes 
with the same n in the coherent state representation ([88l) implies the asymptotic formula for the vertex amplitude to 
be given by the unbalanced square of the above formula [167], namely 
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One can write the previous expression as 
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where 

a<b 

is the Regge like action for Xjjab = Aab the ten triangle areas (according to the LQG area spectrum [l|,f3|). Remarkably, 
the above asymptotic formula is also valid for the case 7 > 1 [8Q]. The first term in the vertex asymptotics is in 
essence the expected one: it is the analog of the 6j symbol asymptotics in three dimensional spin foams. Due to their 
explicit dependence on the Immirzi parameter, the last two terms are somewhat strange from the point of view of 
the continuum fi eld t heoretical view point. However, this seems to be a peculiarity of the Riemannian theory alone 
as the results of [l66j | for the Lorentzian models show. Non geometric configurations are exponentially surpressed 



4. Lorentzian EPRL model 



To each solution one can associate a second solution corresponding to a parity related 4-simplex and, consequently, 
the asymptotic formula has two terms. It is given, up to a global sign, by the expression 
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wher e are constants that do not scale. Non geometric configurations are exponentially surpressed 

In 1713 Freidel and Conrady gave a detailed description of the coherent state representation of the various spin 
foam models described so far. In particular they provided the definition of the effective discrete actions associated 
to each case whic h we presented in (jl01[) . This provides the basic elements for setting up the asymptotic analysis 
presented in (l70j | (the first results on the semiclassical limit of the new spin foam models) which is similar to the 
studies of the asymptotic of the vertex amplitude reviewed above but more general in the sense that the semiclassical 
limit of a full spin foam configuration (involving many vertices) is studied. The result is technically more complex 
as one studies now critical points of the action associated to a coloured complex which in addition of depending on 
group variables g it depends on the coherent state parameters n. The authors of [l70j | write Equation (|10ip in the 
following way 
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They show that those solutions of the equations of motion of the effective discrete action that are non geometric (i.e. 
the contrary of Regge like) are not critical and hence exponentially suppressed in the scaling jf — ^ Xjf with A >> 1. 
If configurations are geometric (i.e. Regge like) one has two kind of contributions to the amplitude assymptotics: 
those coming from degenerate and non-degenerate configurations. If one (by hand) restricts to the non-degenerate 
configurations then one has 
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where rig, and Uf denote the number of edges, vertices, and faces in the two complex A* respectively. There are 
recent works by M. Han w here asymptotics of general simplicial geometry amplitudes are studied in the context of 
the EPRL model [l73 . [T75j . 

The problem of computing the two point function and higher correlation functions in the context of spin foam has 
received lots of attention recently. The framework for the definition of the correlation functions in the background 
independent setting has been generally discussed by Rovelli in [176] and correspods to a special appli cation of a more 
general proposal investigated by Oeckl [I77l4l84 |. It was then applied to the Barrett-Crane model in jl85l4l87 |, where 
it was discovered that certain components of the two point function could not yield the expected result compatible with 
Regge gravity in the semiclassical limit. This was used as the main motivation of the weakening of the imposition 
of the Plebanski constraints leading to the new models. Soon thereafter it was argued that the difficulties of the 
Barrett-Crane model where indeed absent in the EPRL model [l88| |. The two point function for the EPRL model 
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was calculated in [18911 ^-^id it was shown to produce a result in agreement with that of Regge calculus [190|, Il9l| in 
the limit 7—^0. 

The fact that, for the new model, the double scaling limit 7 — > and j — > 00 with 7j=constant defines the 
appropriate regime where the fluctuation behave as in Regge gravity (in the leading order) has been further clarified in 
[192] . This indicates that the quantum fluctuations in the new models are more general than simply metric fluctuations. 
The fact the the new models are not metric at all scales should not be surprising as we know that the Plebanski 
constraints that produce metric general relativity out of BF theory has been implemented only semiclassically (in the 
large spin limit). At the deep Planckian regime fluctuations are more general than metric. However, it not clear at 
this stage why this is controlled by the Immirzi parameter. 

All the previous calculations invol ve a complex with a single four-simplex. The first computation involving more 
than one simplex was performed in |l87l Il93j | for the case of the Barrett-Crane model. Certain peculiar properties 
were found and it is not clear at t his s tage whether these issues remain in the EPRL model. Higher order correlation 
functions have been computed in |l94| . the results are in agreement with Regge gravity in the 7^-0 limit. 
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